A quantum phase transition from triangular to stripe charge order in NbSe2 
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The competition between proximate electronic phases produces a complex phenomenology in 
strongly correlated systems. In particular, fluctuations associated with periodic charge or spin mod- 
ulations, known as density waves, may lead to exotic superconductivity in several correlated ma- 
terials. However, density waves have been difficult to isolate in the presence of chemical disorder, 
and the suspected causal link between competing density wave orders and high temperature super- 
conductivity is not understood. Here we use scanning tunneling microscopy to image a previously 
unknown unidirectional (stripe) charge density wave (CDW) smoothly interfacing with the famil- 
iar tri-directional (triangular) CDW on the surface of the stoichiometric superconductor NbSe2. Our 
low temperature measurements rule out thermal fluctuations, and point to local strain as the tuning 
parameter for this quantum phase transition. We use this discovery to resolve two longstanding de- 
bates about the anomalous spectroscopic gap and the role of Fermi surface nesting in the CDW phase 
of NbSe2. Our results highlight the importance of local strain in governing phase transitions and 
competing phenomena, and suggest a new direction of inquiry for resolving similarly longstanding 
debates in cup rate superconductors and other strongly correlated materials. 



A. INTRODUCTION 

While a classical phase transition separates two states 
of matter at different temperatures, two ordered ground 
states of a material at zero temperature are sepa- 
rated by a quantum critical point (QCP). The competi- 
tion between proximate ordered phases near the QCP 
can dramatically influence a large region of the phase 
diagram [1 J. While the fluctuations from competing 
quantum states lead to exotic physics even at higher 
temperatures, low temperature studies of these states 
can lead to a better understanding of the root of the com- 
petition. Density waves - charge or spin ordered states 
of collective origin driven by instabilities of the Fermi 
surface (FS) - exist in close proximity to superconductiv- 
ity (SC) in several classes of correlated materials fZ*-^, 
and various proposals have recently emerged to study 
their interplay in the presence of strong inhomogene- 
ity in these systems[5]. In this light, it is surprising that 
charge density waves (CDWs) are not fully understood 
even in the weakly correlated and stoichiometric tran- 
sition metal dichalcogenides (TMDCs). While a clas- 
sic CDW arises from strong FS nesting, resulting in a 
sharply peaked susceptibility and a Kohn anomaly at 



the CDW wavevector, the quasi-2D TMDCs are known 
to deviate from this picture|6|. 

2i^-NbSe2 is a layered TMDC which has generated 
much recent interestlZHgl as a model system for under- 
standing the interplay of the CDW and SC phases which 
onset at Tcdw 33 K and Tsc ^ 7.2 K respectively (lOl 
11 J. Despite extensive studyll8l [T2HT5l , several key facts 
about its familiar tri-directional (3(5) CDW remain un- 
resolved, including the role of FS nesting in determin- 
ing its wavevector qsq, and the magnitude of the spec- 
tral gap and its role in the energetics of the transi- 
tion. First, angle-resolved photoemission (ARPES) stud- 
ies have been unable to uniquely identify gag -nested FS 
regions in 2i^-NbSe2||8| [El HSflH. Meanwhile, recent 
studies indicate a broadly peaked susceptibility and a 
soft phonon over a range of wavevectors around qsQ (SI 
[141 [151 , suggesting instead that the g-dependence of 
the electron-phonon coupling could play an important 
role in driving the transition. Second, kinks in tunnel- 
ing spectra at ±35 mV (isx)/ historically identified as 
gap edges, correspond to an anomalously large energy 
scale for the corresponding Tcdw (2£k/3.5/cbTcdw ^ 
7. 05) 1 20 1, while recent ARPES studies indicate a much 
smaller - 3 - 5 mV gap|[T2l[T5l. 
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Our discovery, by low-temperature scanning tunnel- 
ing microscopy (STM), of a IQ CDW with distinct wave- 
length and tunneling spectra from the 3Q CDW, in con- 
junction with band structure calculations, allows us to 
resolve both longstanding questions of the wavevector 
and the gap. First, the distinct wavelengths demonstrate 
that FS nesting plays a negligible role in setting their 
magnitude. Second, the distinct tunneling spectra of the 
IQ CDW region help us disentangle the 3Q CDW spec- 
tra to expose a particle-hole asymmetric gap, riding on 
top of a strong inelastic background. 



B. RESULTS 



Fig. 14 shows a topographic image of a locally com- 
mensurate (3ao) CDW. Its microscopic 3Q nature is con- 
firmed by the existence of a secondary CDW peak in 
the Fourier Transform (FT) in Fig. Ip , in contrast to 
bulk measurements! 9 1. Phase slips result in an overall 
periodicity of Xsq ^ 3.05 ao, corresponding to qsq ^ 
0.328 Qo, where Qo is the Bragg vectorEl [lOl |^_Our 
primary experimental discovery is shown in [Fig. 2 4, 
where regions of unidirectional (IQ) CDW with unique 
wavevector qiq along a single 3Q direction form an 
atomically smooth interface with the 3Q CDW. The ab- 
sence of atomic lattice discontinuities rules out the pos- 
sibility of a NbSe2 poly type interface ||22|. While other 
TMDCs are known to exhibit a thermally induced tri- 
clinic CDW state that varies with doping near Tcdw12TI 
1231 , no such anisotropy has been reported in 2i^-NbSe2. 
Moreover, our observations are at temperatures T <C 
TcBWr where thermal fluctuations are insignificant, im- 
plying that the IQ CDW is a distinct quantum phase. 

The IQ CDW regions take the form of elevated to- 
pographic ribbons, exemplified in [Fig. 2 B-C, suggesting 
a strain-induced phase ( |§ SI The observation of Y- 
junctions between ribbons with differently oriented qiQ 
rules out extrinsic uniaxial strain and suggests instead 
locally varying strain, perhaps due to underlying de- 
fects causing nanoscale buckling of the top few atomic 
layers. From a survey of several ribbons, we place up- 
per bounds of 3% on the vertical strain and 0.06% on the 
lateral strain (see § SI 



Fig. 3A-B show the dominant Fourier peak for the IQ 
modulation, qiQ 2/7 Qo/ corresponding to a wave- 
length, AiQ 0^ 3.5 ao- No similar periodicity or rich 
harmonic structure has thus far been reported in any 
TMDC system[ll]. We develop a phenomenological 
understanding of the IQ harmonic structure following 



McMillan's Landau theory ||2T1 |24l. Rather than a uni- 
form 3.5 ao charge modulation, the system could lower 
its energy by locking the charge modulation to the lat- 
tice with 3 ao periodicity. This would require compensa- 
tion by a one atom phase slip every two oscillations, cor- 
responding to a 27r/3 discommensuration, as shown in 



Fig. 32[21J. The resulting harmonic structure shown in 
Fig. 3D reproduces all observed peak positions. More- 
over, the rich harmonic content we observe is another 
indication of the strong coupling of the electronic mod- 
ulation to the lattice. An even better agreement with 
relative peak heights could be obtained by considering 
spatial variations in the order parameter amplitude [21 J. 

The stark contrast between our observation of two 
CDW wavevectors qiQ and qsq of same orientation but 
13% difference in magnitude, and the recent X-ray mea- 
surements reported by Feng et al.[9\, provides strong 
evidence against FS nesting at one particular wavevec- 
tor as a driving force for either CDW. While our 13% 
wavevector difference arose from moderate anisotropic 
strain (up to 0.06% in-plane), Feng et al. applied hydro- 
static pressure sufficient to induce in-plane lattice distor- 
tions up to 1.6%, yet observed no measurable deviation 
of the CDW wavevector from ^3q|9|. The observed in- 
sensitivity of qsq to hydrostatic pressure would clearly 
indicate that the FS does not qualitatively change in the 
presence of even relatively large lattice distortions, and 
would thus rule out a change in the FS as the source 
of our observed 13% wavevector difference. Further- 
more, consistent with our experiment and with previ- 
ous calculations!^, ^5 1, we find no sharp peak in the 
susceptibi lity (|Fig. 4[ ) computed from our modeled band 
structure ( |§ SI Therefore, our observations and cal- 
culations both indicate that the FS can play only a minor 
role in determining CDW wavevectors in NbSe2. This 
highlights the key role that the g-dependence of alterna- 
tive mechanisms such as electron-phonon coupling may 
play in driving the transition, and particularly the man- 
ner in which these mechanisms may be influenced by 
local strain. 



The two CDW regions display quite different tunnel- 
ing spectra, as shown in [Fig. 5A -B. We employ a fit 
to the NbSe2 band structure ( § SI III) , and impose a 
CDW wavevector q, gap A, and broadening parame- 
ter f, to calculate the density of states (DOS), in Fig. 5 2 
( |§ SHY\ . For the 3Q CDW state the measured dl/dV 
spectrum, proportional to the DOS, is best reproduced 
using q = (0.333 ± 0.004) Qo (the observed local CDW 
periodicity), A = (12 ± 2) mV (which has not been pre- 
viously apparent from direct observations by spectro- 
scopic techniques), and f = 5 mV. The calculations cap- 
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Figure 1: STM of the 3Q CDW. (A) Atomic resolution topograph showing the ^ 3ao periodic CDW. Yellow lines are overlaid on 
a one atom shift of the CDW maximum (phase slip). Inset shows the crystal structure of 2if-NbSe2, with alternating layers of Nb 
(green) and Se (orange) atoms. (B) Fourier Transform (FT) of a larger (^ 45 nm) topograph, displaying a 3Q CDW. Primary CDW 
wavevectors, qi and q2 (yellow circles), a secondary CDW wavevector, ql + q2 (blue circle) and Bragg vectors (dashed red circles) 
are indicated. (C) Average dl/dV spectrum acquired from the area in A, showing kinks at ^ ±35 mV (sk), the minimum Smin 
offset from (Vbias = 0)by^ 17 mV, and marked asymmetry about sf- Setpoint parameters: sample bias, Vsampie = — 50mV(A), 
-60 mV (B), -80 mV (C); junction resistance, i?j - 2.5 GQ (A), 0.1 GQ (B), 0.4 GQ (C); and RMS lock-in excitation, V^od - 3.5 mV 
(C). 




Figure 2: The IQ — 3Q Interface. (A) Topograph showing an atomically smooth interface between the 3Q (bottom) and IQ (top) 
CDWs. (B) Topograph showing formation of IQ CDW along a 'ribbon'. Ribbons are typically 10 — 20 nm in width, and elevated by 
20 — 40 pm. Dashed yellow lines indicate the approximate extent of the ribbon; the red arrow points along the IQ CDW direction 
(Qiq)- (O Topograph showing a F-junction between three such IQ ribbons, each with a different IQ CDW direction. The 
topograph in A has been leveled by removing a polynomial background to clearly show the CDW interface. Setpoint parameters: 
^sample = "50 uiY ) Rj = 1 GQ (A,C), 5 GQ (B). 



ture the overall shape, width, and center Smin of the gap 
structure within 30 mV of ( |§ SI IV\ . The fact that ^min 
is offset from should be unsurprising for a quasi-2D 
system 1251 , but had not been understood or observed 
in NbSe2 until now, due to limitations of spectroscopic 
techniques which are sensitive to filled states only 1121. 

We disentangle the CDW gap from other effects in the 
3Q spectra through a comparison with the IQ spectra 
in |Fig. 5p . These F-shaped IQ spectra resemble the lin- 
ear tunneling conductance background historically at- 
tributed to the inelastic coupling of tunneling electrons 
to a flat bosonic spectrum ! 26 1 . That this background 



is much stronger in the IQ region, obscuring band 
structure effects, is likely a strain-induced phenomenon, 
which may be related to the buckling and associated de- 
coupling of the topmost layers in the IQ region. Mean- 
while, present in both IQ and 3Q spectra (thus unlikely 
to be associated with these different CDWs), yet absent 
in calculations (thus unlikely to be a band structure ef- 
fect), are the ±35 mV kinks, previously and mistakenly 
identified as the CDW gap|20|. We universally observe 
the kinks even well above Tcdw/ which further demon- 
strates their lack of bearing on the CDW phase ( |§ SI TV) . 
ARPES studies observe a prominent band structure kink 
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Figure 3: The IQ CDW. (A) FT of the out-of-feedback current at +50 mV over a IQ region (upper left quadrant of [Fig. 2 4). The 
dominant CDW wavevector (yellow circle) and Bragg vector (dashed red circle) are indicated. (B) Linecut of the FT intensity 
parallel to the red line in A from the center to the Bragg peak, in units of the Bragg vector Qq. The dominant peak, qiq ^ 2/7Qo 
(solid blue arrow) and its harmonics (dashed blue arrows) are identified, and are distinct from the 3Q wavevector qsq (dashed 
green line). The remaining peaks are Bragg reflections of these three peaks. Setpoint parameters: V^sampie = —50 mV; Rj — 0.2GQ. 

(C) Phenomenological model of the observed IQ CDW wavelength, Aiq :^ 3.5 ao (details in text). The atomic periodicity is 
indicated by green circles and the phase of the CDW order parameter from the model is shown in blue[21J. A topographic linecut 
(red dots) is extracted from Fig. 2j 4 along the red arrow {qiq), filtered to remove atomic corrugations, and overlaid for comparison. 

(D) Simulated FT intensity from the cartoon CDW modulation in C, for comparison with experimental peak positions in B. 
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Figure 4: Susceptibility. The non-interacting susceptibility 
Xo(^, 00 — 0) calculated from a fit to the 2if-NbSe2 band struc- 
ture \% SI III^ along the V — M direction of the reciprocal lat- 
tice, displaying a broad maximum over a range of wavevec- 
tors: q - 0.25 - 0.4(6l[ll[ia|. The CDW wavevectors qiq, q^q 
and (5o/3 are overlaid for comparison. Inset shows the Bril- 
louin Zone (BZ) of the hexagonal 2i7-NbSe2 lattice. 

at a similar energy in the Se r-pocket|[l5l[l8l, attributed 
to coupling to an optical phonon ||T3l . We therefore con- 
clude that this self-energy effect is responsible for the syl 
kinks in the tunneling spectra as well. The discrepancy 



between the data and band structure calculations above 
~ 30 mV in |Fig. 5p can thus be attributed to the inelastic 
tunneling background and self -energy effects. 



C. DISCUSSIONS 

We therefore resolve a longstanding debate about 
the anomalous CDW gap magnitude reported by 
STM measurements[20J, and caution that not all e^- 
symmetric kinks in tunneling spectra are associated 
with order (e.g. density wave or superconducting gaps). 
On the contrary, we emphasize that the true CDW 
signature in NbSe2 is offset from sy, which has con- 
fused an active research community for two decades, 
and has been disentangled now only by a combina- 
tion of spatially resolved filled and empty state spec- 
troscopy of a proximate (IQ) phase, and band struc- 
ture calculations 1 25 1 . This emphasizes the need for full 
experimental exploration of proximate phases in other 
pertinent materials, combined with quantitative mod- 
eling. We further suggest that controlled local strain, 
through epitaxy or intentional defects, may be a use- 




Figure 5: Spectroscopy across the IQ 

red arrow in 



3Q Interface. (A) Linecut of dl/dV spectra across the IQ — 3Q interface taken along the 



Fig. 2 4. (B) Representative spectra in the 3Q (red) and IQ (green) regions of 



Fig. 2 4 overlaid for comparison. The 



IQ spectra have a minimum close to £f, a deep V-shape with reduced asymmetry, and kinks at ±£k- The spectra are normalized 
at —50 mV. Setpoint parameters: Vsampie = —50 mV; Rj — 0.2 Gil; Vmod = 3.5 mV. (C) Calculated DOS for NbSe2 in the 'normal' 
state (black) and in the presence of a 3Q CDW (blue) using q = 0.333 Qo, A = 12 mV, f = 5 mV, compared with the 3(9 STM 
spectrum (red) (details in§SI IV). The calculations reproduce the observed asymmetry, offset Smin, and shape of the gap structure. 



ful tuning parameter to access the necessary proximate 
phases for comparison. 

Beyond providing new^ insight into the nature of the 
3(5 CDW in NbSe2, our v^ork motivates the utility of the 
IQ — 3(5 interface in NbSe2 as a platform to explore com- 
peting quantum phases in the weakly correlated limit, 
as a step towards understanding them in strongly cor- 
related systems. In the Landau picture of CDWs[24J, a 
quantum phase transition between 3(5 and IQ states can 
arise by tuning the coefficient of the interaction term be- 
tween the three inequivalent CDW propagation direc- 
tions (though in our case, the magnitude of q differs 
between the two states). In NbSe2, even at low tem- 
peratures T <C TcDW/ where the amplitude of the or- 
der parameter is already large, moderate strain is seen 
to have a strong influence, indicating that the system 
is intrinsically close to the QCP separating these states. 
We note that a related phase transition between the ob- 
served IQ CDW phase and a 'hidden' 2Q phase has been 
suggested, but not directly visualized, in the rare-earth 
tritelluridesl27n28J. 

Our discovery provides a new perspective on the role 
of density wave order in complex systems. First, our res- 



olution of two longstanding debates about NbSe2 puts 
this much-studied material on firmer footing as a well- 
understood model system for CDW studies and compet- 
ing ground states in superconductors. We have disen- 
tangled the true CDW gap, and clarified that FS nesting 
plays a minor role in determining the CDW wavevec- 
tors in this material, thereby highlighting the role of 
other mechanisms in driving the transition. Second, our 
revelation of a particle-hole asymmetric CDW gap em- 
phasizes the limitations of filled-state-only probes, e.g. 
ARPES, for investigating phases other than SC - which 
is unambiguously particle-hole symmetric. Full spectral 
probes such as STM, in combination with quantitative 
calculations, are necessary to understand the competi- 
tion between SC and particle-hole asymmetric phases. 
Third, our observation of the local effect of even mod- 
erate strain in driving a quantum phase transition calls 
for a reinvestigation of possible phase inhomogeneity 
in other strongly correlated systems, where larger strain 
may occur 1 29 . 30] . 

In the cuprate superconductors, an analogous 
phase boundary between unidirectional (1(5) charge 
'stripes' and bidirectional (2(5) 'checkerboard' has been 
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predicted ||3T}|32|. The introduction of quenched disor- 
der results in discommensurations in the 2Q phase and 
disordered orientational order in the IQ phase, mak- 
ing them hard to distinguish - especially in the cuprate 
BSCCO, thought to be in proximity to the IQ — 2Q 
phase boundary 1311 l32l . Recent STM studies of the 
~ 4 ao charge order in BSCCO have had conflicting in- 
terpretations, with independent suggestions of fluctu- 
ating 2Q and IQ order|33. 34 1. However, the influ- 
ence of strain, from the supermodulation lattice buck- 
ling, or from randomly distributed dopants, is seldom 
accounted for. Previous studies have shown that both 
these strain phenomena correlate with nanoscale elec- 
tronic inhomogeneity 1291 l35l . A possible explanation is 
local stabilization of the IQ state, producing 1Q—2Q and 
IQ — IQ interfaces with spectral differences, analogous 
to Fig. 2 and Fig. 5A-B. While the presence of strong dis- 



order (up to 12% strain variations on a nanometer length 
scale 1 29 1) complicates the interpretations in BSCCO, we 
stress the importance of isolating and modeling strain 
effects for better understanding and control of the phase 
transitions in cuprates. Finally, the microscopic visu- 
alization of the role of strain in stabilizing new order 
suggests a controlled route towards engineering novel 
quantum phases and interfaces and studying symmetry 
breaking in strongly correlated materials. In this regard, 
we suggest a connection to the emerging importance of 
strain as a route to high-Tc superconductivity in novel 
iron-based materials I l36l |37| . 



D. METHODS 



STM Experiments. Measurements were performed 
using a home-built STM at temperatures between 2 — 10 
K. Magnetic fields of up to 5 T were used to suppress 
the superconducting state as needed. Single crystals of 
2i^-NbSe2 were cleaved in situ in cryogenic ultrahigh- 
vacuum and inserted into the STM. A mechanically 
cut Ptir tip, cleaned by field emission and character- 
ized on gold, was used for the measurements. Spec- 
troscopy data were acquired using a lock-in technique 
at 1.115 kHz. The topographic and spectroscopic signa- 
tures of the IQ ribbons have been verified with several 
tips. 

Band Structure and DOS Calculations. The band 
structure of 2i^-NbSe2 close to consists of two Nb—4d 
derived bands and one Se—4p 'pancake'-shaped hole 
pocket! 6 1. The Nb— 4(i bands are modeled using a tight- 
binding fit to the observed ARPES band structure] 15], 
and the Se— 4p pocket is approximated by a simple 



quadratic form to fit LDA calculations f6l (details in 
|§ SI III) . Using all three bands, the DOS in the presence 
of a CDW is calculated by imposing a constant coupling 
between electronic states connected by any one of the 
three g'-vectors. The strength of the coupling A is taken 
as a free parameter in the reproduction of the experi- 
mentally observed DOS, and the size of q is allowed to 
vary slightly around the observed value of 0.328 Qo (de- 
tails in |§ SI IVp . 
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SI I. SAMPLE CHARACTERIZATION 
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Figure SI: Sample Characterization. (A) The layered hexag- 
onal crystal structure of 2if-NbSe2, with alternating sand- 
wiches of Se-Nb-Se. (B) Temperature dependence of the re- 
sistance for the sample batch used for this study, showing a 
superconducting transition at 7 K. 



2i^-NbSe2 is a layered transition metal dichalco- 
genide with a hexagonal structure and Dq^ space group 
symmetry The unit cell ( |Fig. SH 4) consists of two sand- 
wiches of Se-Nb-Se. The crystal typically cleaves be- 
tween the neighboring Se layers, coupled by weak van 
der Waals forces. 

Single crystals of 2iJ-NbSe2 were grown by chemical 
vapor transport using iodine as the transport agent. A 
transport characterization of the sample batch used in 
this work is shown in |Fig. Sip . The superconducting 
transition is observed at Tsc ~ 7K. The residual re- 
sistivity ratio (RRR), defined as the ratio of resistances 
i?(295K)/i?(7.5K),is- 16. 



SI II. IQ 'RIBBONS': HEIGHT AND ORIENTATION 



The IQ CDW typically appears in regions which per- 
sist in one direction with apparent 20 — 40 pm topo- 
graphic elevation, forming a 10 — 20 nm wide 'ribbon' 
structure ( |Fig. S2| 4-C). We note that the topographic el- 
evation ^stm(^, ^0, ^o) as measured by maintaining a 
constant current Iq with bias setpoint Vq at lateral tip 
position r = {x^y) can be given by 



Here zrir) is the true topographic corrugation of the 
sample, hi{f) is a measure of the local tunnel barrier 
height (LBH), and D (r, eV) is the local density of states 
(LDOS) of the sample at energy eV. Because of the log- 
arithmic sensitivity of ^stm(^, ^o, ^o) to the integral of 
the LDOS from the Fermi energy, (corresponding to 
V = 0)up to the bias setpoint eVor STM topographs may 
contain electronic artifacts masquerading as geometric 
effects. Therefore, we present two pieces of evidence for 
the true geometric elevation of these ribbons. 

First, a tabulation of the relative orientation Or^q^ of 
the ribbon to the nearest Bragg vector of the underlying 
hexagonal lattice for the various ribbons imaged in the 
study shows a seemingly random spread from —30° to 
30° - the full range of available angles ( |Fig. S2p -£). Fur- 
thermore, these ribbon structures can intersect to form 
X as well as Y junctions ( |Fig. S2 4-.B), and the angle be- 
tween intersecting ribbons varies from 40° to 60°. The 
fact that ribbon orientation does not respect lattice sym- 
metry strongly suggests a true geometric, rather than 
electronic origin of their apparent height. We contrast 
this observation with enhanced STM topographic cor- 
rugation associated with predominantly electronic fea- 
tures in a wide range of other materials, which respect 
the symmetry of the hexa gonal 1381 l39l or square L33 J lat- 
tice. 

Second, the measured height of these ribbons in STM 
topographs exhibits < 5% dependence on bias setpoint 
within 400 mV below the Fermi energy ( |Fig. S3 4-C). 
We note that this energy range over which the mea- 
sured height of the ribbons is invariant is much larger 
than the spectral range of CDW variation in the DOS 
50 mV). We further note the contrast between the 
bias-independent ribbons, and single atom impurity res- 
onances, whose measured 'height' varies by 50 — 70% 
between Fig. S3| 4 and B. Therefore we conclude that the 
measured height (20 — 40 pm) and width (10 — 20 nm) of 
these ribbons has a predominantly geometric origin. 

Having established the topographic origin of these 
ribbons ( Fig. S2| and |Fig. 53) , we suggest that these rib- 
bons are likely a topographic rippling of the top few lay- 
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Figure S2: IQ 'Ribbons': Orientation. (A, B) STM topographs showing 'ribbons' of unidirectional (IQ) CDW, intersecting to 
form an X-junction (A) and a F-junction (B), with CDW wavevector qiQ varying between the arms. Dashed yellow lines indicate 
the approximate extent of the ribbons. Setpoint parameters: 14ampie = — 50mV for both; i?j = 5 (A), 1 GQ (B). (C) STM 
topograph of a IQ — 3Q interface (top arm of X-junction, also shown in Fig. 2p ), demonstrating the definitions of the relative 
orientations between the ribbon (green arrow) and the nearest Bragg vector (dashed red arrow), and ^R,qig, between the 

ribbon and IQ wavevector qiq (solid red arrow). (D) A table detailing the values of ^r,Qq and fe,qiQ observed in the various IQ 
ribbons studied in this work, with the first entry corresponding to C. (E) A visual illustration of the spread of values in table D. 
The dark blue (single), green (X-junction) and cyan (F-junction) lines describe the orientation of the various ribbons with respect 
to the IQ wavevector qiQ along the (5o,i Bragg vector (solid red arrow) and another Bragg vector Qo,2 (dashed red arrow). 



ers. These ribbons may arise during the cleaving process 
due to underlying growth defects w^hich can intercalate 
betw^een Se-Nb-Se sandw^ich layers. We note that simi- 
lar topographic ribbon deformations have recently been 
observed in another layered chalcogenide (Bi2Te3)[40j. 

We estimate the in-plane and out-of-plane lattice 
strain associated with the topographic ribbon features. 
Using the maximum topographic elevation of an ob- 
served ribbon (40 pm), we can put an upper limit on 
the out-of-plane distortion by assuming that a mini- 
mum of 2 sandwich layers are elevated (any fewer, and 
the defects causing the elevation would be likely visi- 
ble in our topographs). The out of plane distortion is 
therefore < 3% (40pm/12.54 A) of the unit cell spac- 
ing. To measure the in-plane distortion, we first use 
the Lawler-Fujita algorithm ic , which can determine 
the lateral location of atoms with precision ~ 2% of the 
lattice spacing|42, 431. With this algorithm, we do not 
observe any change in the lattice constant across the rib- 



bon, which places a direct experimental upper limit on 
the in-plane distortion of ~ 2%. 

However, we can estimate the actual in-plane distor- 
tion indirectly from the measured out-of-plane distor- 
tion, by modeling the ribbon as a half-period of a si- 
nusoid with height h (40 pm) and width W (10 nm) 
( |Fig. S3p ). The total lateral deformation due to such a 
ribbon is 5W ~ 45 pm, corresponding to ^ 0.06% of the 
lattice spacing. As this is well below the resolution of 
the Lawler-Fujita algorithm, it is not surprising that the 
in-plane distortion is not detectable in STM topographs. 
From the upper bounds of 3% on the vertical strain and 
0.06% on the lateral strain, we note that the magnitude 
of the strain field leading to the formation of these rib- 
bons is moderate, in comparison to some other corre- 
lated materials lE9l l3Ql . We also note that while the mag- 
nitude of lattice distortion of these ribbons may seem 
small in the context of the observed quantum phase in- 
homogeneity, a comparison with other known materials 
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Figure S3: IQ 'Ribbons': Bias Setpoint Dependence. (A, B) Topographs of IQ ribbons over the same spatial area acquired 
with different bias setpoints: -100~mV (A) and -400 mV (B), with i?j = 10 in both. The dashed yellow lines indicate the 
approximate extent of the ribbons, and the dashed green circles enclose triangular impurities, visible in B with a 50 — 70% larger 
apparent height than in A. (C) Linecuts taken through the topographs in A (red) and B (blue) transverse to the lower ribbon, 
along the dashed white arrows. The measured ribbon height varies less than 5% between bias setpoints -100 mV and -400 mV. 
The linecuts have been laterally averaged over a 3 nm width indicated by the dashed grey lines. (D) A cartoon representation of 
the lattice distortion caused by the formation of such a ribbon, modeled as a half period of a sinusoid. The red circles correspond 
to rows of displaced Se atoms, while the dashed red circles represent their original undistorted positions. The ribbon has height 
h (40 pm) and width W (10 nm), resulting in a total lateral distortion of 5W ^ 45 pm across the ribbon. 



suggests that strain of this magnitude can be sufficient to 
drive a transition to the unidirectional CDW phase 1 28 1. 

We previously discussed the ribbon orientation with 
respect to the lattice (6>r^Qq); we now consider the rib- 
bon orientation with respect to the IQ CDW wavevec- 
tor (^R,giQ), also detailed in Fig. S2p -£. In a simple pic- 
ture of the strained ribbon structure, we would expect 
the ribbon-induced strain to couple strongly to the IQ 
CDW orientation either parallel or perpendicular to the 
ribbon, and thus we would expect to observe values of 
^R,qiQ either between - 30° or between 60 - 90° . Yet we 
often find O^^q^Q to be in the 30 — 60° range as well. In- 
sufficient statistics prevent us from inferring a clear con- 
nection between ribbon orientation and qiq orientation, 
but the wide distribution of relative angles suggests the 
complexity of the interaction. 



(S2) 



With the values of A = -5.4 eV and B = -0.65 eV, 
this model accurately reproduces the dispersion ob- 
tained in LDA calculations by Johannes et al\6\, as 
shown in Fig. S4| 4. To model the Nb-4(i bands, we use 
a tight-binding fit to the band structure observed in 
ARPES experiments, as reported by Rahn et al\15\. We 
find that a small (~ +16meV) offset in the chemical 
potential relative to the parameters used by Rahn et al. 
was needed to reproduce the observed DOS from STM 
measurements. This offset is within the accuracy of the 
tight-binding fitting scheme] 15 1. The band energies of 
the tight-binding description are given by 1151 



E = tQ^ti[2 cos ^ cos T] + cos 2^] 

Sim. BAND STRUCTURE CALCULATIONS 

+ t2 [2 cos 3^ cos 77 + cos 277J 
+ [2 cos 2^ cos 2r] + cos 4^] (S3) 
The band structure of 2i^-NbSe2 close to consists + ^4 [cos ^ cos 3/? + cos 5^ cos r] + cos 4^ cos 2r]\ 

of two Nb-4(i derived bands and one Se-4p 'pancake'- ^ j2 ^qs 3^ cos 3/^ + cos 6^] 

shaped hole pocket||6l |17|. Close to e^, the Se-4p 

pancake-shaped hole pocket which surrounds the F- 

point can be modeled by a simple quadratic form, where ^ = \kxaQ and r] = ^y/^kya^ and ky is along 
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Figure S4: Band Structure Calculations. (A) The calculated band structure of 2if-NbSe2 along high-symmetry directions, show- 
ing the two Nh-M bands obtained using an ARPES tight-binding fit| 15 1, and the Se-4p band, modeled as a parabolic fit to LDA 
calculations of Johannes et aL\6\, compared with the LDA calculations (dashed black lines). (B) The Fermi surface obtained using 
this band structure fit, with the BZ shown in black. 



r — M. The values of the tight binding parameters (in- 
cluding the offset) used in this work are: 



Parameter 


Nb Band (1) 


Nb Band (2) 


to 


26.9 


219.0 


ti 


86.8 


46.0 


t2 


139.9 


257.5 


u 


29.6 


4.4 


tA 


3.5 


-15.0 




3.3 


6.0 



Table SI: Tight-binding parameters for the two Nb bands (in 
meV). 



SI IV. DENSITY OF STATES CALCULATIONS 



To calculate the DOS in the presence of a SQ CDW, we 
impose a coupling between states connected by any one 
of the three ^-vectors, given by: 



^CDw = A.^ (4^,c,- + h.c.) 



(S4) 



The strength A of the coupling is taken as a free param- 
eter in the reproduction of the experimentally observed 
DOS, with a broadening parameter fixed at F = 5 meV. 
Adjusting the size of q slightly around the observed 
value of = 0.328 Qq, we find the best match with 
STM spectra using q = 0.333 Qo (corresponding to the 
locally commensurate CDW periodicity), A = 12 meV. 



With these parameter values, the gap structure in the 
calculated DOS closely approaches the overall shape, 
width and center of the gap structure seen in the mea- 
sured data within ±30 meV of the Fermi energy sy, as 
shown in Fig. 5p . 



To demonstrate the accuracy of this fit, we show the 
effects of varying the wavevector q b y 1 — 2% ( |Fig. S5 4- 
B), and the gap value A by 10 - 20% ( [Fig. S5| 2-D). Using 
these fit parameter variations, we estimate the errors for 
q and A to be 0.004 Qo and 2 meV respectively. The value 
of A may however be an overestimate, leading to a sys- 
tematic error of the same order as the fit uncertainty, 
since the described procedure does not take into account 
the particle-hole symmetric inelastic background in the 
experimental DOS. Accurate modeling of the inelastic 
background would require detailed temperature depen- 
dent spectroscopic data, which is beyond the scope of 
this work. 

Crucially, we note that the particle-hole asymmetry 
in the CDW gap, with its minimum centered above sy, 
cannot be removed by the subtraction of a particle-hole 
symmetric background. Likewise, the striking deviation 
of our fitted gap parameter (A = 12 meV) from previous 
results (four times larger than the 3 meV value detailed 
by Borisenko et al. W2\ . and three times smaller than the 
35 meV value detailed by Hess et a/.f20|) far exceeds fit 
or systematic uncertainties. 

To experimentally verify the lack of bearing of the 
35 mV kinks on the CDW phase, we performed spec- 
troscopy up to 45 K, and universally observed the pres- 
ence of these kinks in the STM spectra well above Tcdw 
(~ 33 K). A comparison of the typical spectrum ac- 
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Figure S5: Variations in DOS Fit Parameters. Calculated DOS spectrum using the band structure fit in the presence of a 3Q 
CDW, showing the effects of varying the fit parameters, wavevector q and gap value A around the best fit values {q = 0.333 Qo, 
A = 12 meV, dark blue), compared with the STM data (red). (A, B) show the effects of varying the wavevector q around 0.333 Qo 
by ±1% and ±2% respectively. (C, D) show the effects of varying the gap value A around 12 meVby ±10% and ±20% respectively. 
Error bars for q (0.004 Qo) and A (2 meV) are deduced using these variations. 
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Figure S6: Spectroscopy above Tcdw- Spectrum acquired at 
40 K (red) compared with that acquired at 2 K, 6 T (blue) on 
the same cleaved surface (within 300 nm). The 2 K spectrum 
has been thermally broadened to 40 K for ease of comparison. 
The spectral kinks at ±35 mV (±£k) are distinguishable well 
above Tcdw, as shown by the guides to the eye (c.f. ¥ig. 1 Z for 
a low temperature comparison). 



quired at 40 K to that acquired at 2 K, is show^n in Fig. S6 



smeared by 0{3kBT), i.e. ~ 10 mV, resulting in a broad- 
ening of the kinks. Despite this, the kinks remain distin- 
guishable, and are present throughout all spatial regions 
studied in this work. 



For completeness, we show^ in Fig. S74 the calcu- 
lated DOS in the presence of a IQ CDW at the exper- 
imentally observed wave vector, qiq = 0.286 Qq, us- 
ing A = 12meV. The lack of correspondence between 



We note that the data acquired at 40 K are thermally 



this calculation (green curve in Fig. S74) and the mea- 
sured dl/dV spectrum (green curve in Fig. 5p ) can be 
attributed to the increased intensity of the inelastic back- 
ground in the buckled region of the IQ ribbons. This 
background is evident in the F-shape of the dl/dV 
spectrum in [Fig. 5p , centered close to sp/ as explained 
theoretically 1 26 1 and observed experimentally across a 
wide variety of materials li26M44rr50 J . 

We also compare our results to a calculation of the 
DOS based on a tight-binding fit to the full three- 
dimensional LDA band structure reported by Johannes 
et alM- The DOS obtained using this 3D LDA fit is 
compared to the results based on the 2D fit in |Fig. S7p . 
We find that the STM data are best reproduced using 
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Figure S7: Comparisons of DOS Calculations. (A) Calculated DOS spectrum using the ARPES tight-binding fit - in the presence 
of a 3Q CDW {q = 0.333 Qo, blue) and a IQ CDW (^iq = 0.286 Qo, green), with A = 12meV. (B) Calculated DOS spectrum 
using the tight-binding fit to ARPES data|15| compared with a tight-binding fit to LDA calculations] 6 1 in the 'normal' state and 
in the presence of a 3Q CDW {q = 0.333 Qo, A = 12 meV). The DOS spectrum calculated using the fit to the ARPES data better 



reproduces the gap feature observed in the STM dl/dV spectrum in Fig. 5 



the band structure observed by ARPES, and that there 
is a noticeable difference between the depths of the 
gaps in the two-dimensional ARPES based and three- 
dimensional LDA-based band structure fits (despite in- 



dependent parameter optimization), which is indica- 
tive of some difference between the surface and bulk 
dispersions (S). 



